Scaling of nano-Schottky-diodes 
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A generally applicable model is presented to describe the potential barrier shape in ultra small 
Schottky diodes. It is shown that for diodes smaller than a characteristic length l c (associated 
with the semiconductor doping level) the conventional description no longer holds. For such small 
diodes the Schottky barrier thickness decreases with decreasing diode size. As a consequence, the 
resistance of the diode is strongly reduced, due to enhanced tunneling. Without the necessity of 
assuming a reduced (non-bulk) Schottky barrier height, this effect provides an explanation for several 
experimental observations of enhanced conduction in small Schottky diodes. 



The effect of downscaling the dimensions of a device on 
its electrical transport properties is an important topic 
today. Extremely small diodes have been experimentally 
realized and characterized in various systems, e.g. car- 
bon nanotube heterojunctions M, junctions between p- 
type and n-type Si nanowires g| or junctions between 
the metallic tip of a scanning tunneling microscope and a 
semiconductor surface || ||. These experiments showed 
several deviations from conventional diode behavior. De- 
spite some modelling in truly one-dimensional systems 
H [| , little work has been done on modelling the effects 
of downscaling a conventional diode, in the regime where 
quantum confinement does not play a role. 

In this paper we present a simple model (based on the 
Poisson equation) describing the barrier shape in a diode, 
that is readily applicable to arbitrarily shaped small junc- 
tions. It is related to descriptions of inhomogeneities in 
the Schottky barrier height (SBH) in large diodes |7j], bar- 
rier shapes in small semiconducting grains [0 and charge 
transfer to supported metal particles |J. Although we 
restrict ourselves to metal-semiconductor junctions, the 
model can easily be adapted for e.g. p-n-junctions. The 
main result is that if the size of the metal-semiconductor 
interface is smaller than a characteristic length l c , the 
thickness of the barrier is no longer determined by the 
doping level or the free carrier concentration, but instead 
by the size and shape of the diode. The resulting thin 
barrier in small diodes will give rise to enhanced tunnel- 
ing, qualitatively explaining measurements of enhanced 
conduction ||, |[ [hJ , without the necessity of assuming 
a reduced SBH. Moreover, experimentally observed scal- 
ing behavior and deviating TV-curve shapes can be 
explained. 

The transport properties of a Schottky diode are gov- 
erned by the potential landscape which has to be tra- 
versed by the charge carriers. First, we study an easily 
scalable and highly symmetrical model system, namely 
a metallic sphere embedded in semiconductor (see Fig- 
ure 1, upper left inset). The radius a of the metallic 
sphere is a measure for the interface size: for large a, we 
expect to find the well known results for a conventional 
diode, while decreasing a gives the opportunity to study 
finite size effects. 



We only model the barrier shape in the semiconductor; 
the SBH ifB is accounted for in boundary conditions and 
is considered as a given quantity. For simplicity, the de- 
pletion approximation [p"l| is adopted, which is valid for 
a wide range of realistic parameters. Moreover, the space 
charge region is assumed to be homogeneously charged, 
an assumption that will be discussed later. Solving the 
Poisson equation in n-type silicon with the boundary con- 
dition that the charge on the sphere cancels the total 
charge in the space charge region, we find for < x < w 
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where x is the radial distance from the interface, w the 
depletion width and L D = y / e s kT/(e 2 Nd) the Debye 
length. The zero-point of the potential is chosen in the 
semiconductor bulk. The value of w is fixed by the sec- 
ond boundary condition V(Q) — V s , where V s is the total 
potential drop over the space charge region and satisfies 
V s = fa - tp B )/e - V (with <p s = E c - Ei). Eq. (jjj) is 
valid for small bias voltage V. The limited validity of 
the depletion approximation at finite temperatures only 
affects the tail of the barrier (where V(x)| < kT), which 
is unimportant for the transport properties. From the 
equation, it can be seen that the characteristic length 
scale of this system is 



l c d = L D ^2eV s /kT 



2e s V s 
eN d 



By comparing the diode size a to l c we can decide whether 
the diode is 'small' or 'large'. In the lower right inset of 
Figure 1 the value of l c is plotted versus doping concen- 
tration 7V(j. 

An important quantity for electrical transport is the 
Schottky barrier thickness. In Figure 1, the barrier 
full width at half maximum (FWHM, X1/2) calculated 
from Eq. (|l|) is plotted as a function of diode size a. 
From the figure it is clear that for a ^> l c the value 
of X1/2 approaches a constant, which was expected for 
a large diode. Indeed, for a » l c , Eq. (lij) reduces to 
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FIG. 1: Plot of the calculated barrier FWHM xi/ 2 as a func- 
tion of diode size a (based on Eq. (|l|)), both in units of l c . The 
dashed lines represent the asymptotic values for a^> l c (con- 
ventional diode) and a <C l c (new regime) respectively. The 
lower right inset is a plot of l c as a function of doping level 
N d in silicon (e s = 11.7) for ip B = 0.67 eV and T = 300 K. 
The upper left inset schematically shows the model system, a 
metallic sphere embedded in semiconductor. 
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FIG. 2: The solid lines are contours of the barrier FWHM 
for various disc-shaped contacts (see inset; radii ranging from 
30 nm (a) to infinite (e)), taken from a numerical solution of 
the Poisson equation in silicon. It clearly shows the contact 
size dependence for contact radii smaller than l c ~ 750 nm. 
The dashed lines are the FWHM-contours of the barrier for 
the three smallest diodes, neglecting the screening effect of 
the semiconductor space charge region. The inset indicates 
the plane of cross-section shown in the figure. 



V(x) = —^^{x — w) 2 , which is the well-known textbook 
PI result for band bending in the depletion approxi- 
mation for an infinitely large diode. Both the depletion 
width w = yj {2e s / eNa)V & and X\n are in that regime 
independent of a. 

Figure 1 shows that for a < l c the value of X\/2 is no 
longer constant, but decreases with decreasing a. For 
a <C l c it approaches X1/2 — a, i.e. the barrier thickness 
equals the diode size. This also follows from Eq. (0), 
which reduces to V{x) = V s ■ a/(a + x) for a <C l c and 
x <C w (that is, close to the interface). Note that this 
is exactly the potential due to the charged sphere only. 
In this regime, the effect of the semiconductor space 
charge on the barrier shape and thickness can be ne- 
glected. This can be understood from the fact that the 
screening due to the space charge region takes place on 
a length scale l c , as in conventional (large) diodes. How- 
ever, from Gauss's law it follows that any charged object 
of typical size d < oo in a dielectric medium gives rise 
to a potential that behaves roughly as V(r) oc d/r. This 
Coulomb potential can be further screened by the for- 
mation of a space charge layer of opposite sign, but that 
additional screening can be neglected if d <C l c . This ob- 
servation holds for any interface with typical dimensions 
much smaller than l c . 

In a geometry that can actually be fabricated, the Pois- 
son equation must be solved numerically. We have done 
this for n-doped silicon (jV<j = 10 15 cm~ 3 ) in contact with 
metallic circular disks of various radii. In all further 
calculations ips — 0.67 eV was used, which is the bar- 
rier height of the CoSi 2 /Si(lll)-interface Q. Figure 2 
shows the FWHM-contours of the barriers as resulting 



from these calculations. Also shown are the FWHM- 
contours of the barrier due to the metallic contacts only, 
illustrating the negligible effect of the space charge region 
on the barrier thickness in very small diodes [p"3| . 

To study the effect of the reduced barrier width on the 
transport properties of a small Schottky diode, a trans- 
mission coefficient T(E, V) was obtained for the barrier 
shape from Eq. (El). This was done in a one-dimensional 
fully quantum mechanical calculation H|. Note that 
T(E, V) is implicitly dependent on temperature and dop- 
ing level, because these quantities influence the position 
of the Fermi-level in the bulk semiconductor. The current 
density is then given by 

/•oo 

J(V) oc / T(E, V) [f(cp s + E)- f( Vs + E + V)} dE, 
Jo 

from which it follows that the zero bias differential con- 
ductance satisfies 



dJ 
dV 



T(E,V)f'(<p.+E) dE. 



Here, / is the Fermi-Dirac distribution function and E 
the energy above the semiconductor conduction band 
edge. Transport due to electrons at energies below the 
barrier maximum (E < V s ) is regarded as tunneling, 
while for E > V s we speak of thermionic emission. Obvi- 
ously, the contribution of thermionic emission is almost 
independent of the barrier thickness, while tunneling is 
strongly dependent on the barrier thickness. 

In Figure 3, the calculated zero bias differential con- 
duction is plotted as a function of diode size a for several 
values of N^- For a > l c this quantity is independent 
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FIG. 3: The contributions of tunneling and thermionic emis- 
sion to the zero bias differential conductance, plotted as a 
function of diode size a for various doping concentrations. 
The vertical arrows indicate the values of l c . The parameters 
are the same as in Figure 2. 




FIG. 4: Calculated IV-curves for various diode sizes. The 
large diode curve has the expected exponential shape. The 
qualitative appearance of the curves changes drastically with 
decreasing diode size. The curves of the larger diodes have 
been scaled vertically. 



of a. For smaller values, the tunnel current starts to in- 
crease rapidly, eventually leading to a strong increase of 
the total conduction. 

Moreover, the shape of the TV-curves changes with de- 
creasing diode size. Our calculations (Fig. 4) show that 
for large diodes the TV-curve has exactly its expected 
exponential shape (7 ex [exp(eV / kT) — 1]). Apart from 
the total current increase, in small diodes the relative 
contribution of the reverse current starts to increase and 
eventually — in extremely small diodes — the reverse cur- 
rent exceeds the forward current, thus reversing the rec- 
tifying behavior of the diode. 

Note that the image charge effect [O was neglected so 



far. However, inclusion of this effect would only enhance 
the phenomenon mentioned above, as it reduces the ef- 
fective barrier height and width even further, especially 
in barriers which are narrow already. 

One more issue that needs to be discussed is that of 
discrete random dopants. In our analysis, the dopants 
played a role in determining the Fermi-level position in 
the semiconductor bulk and were considered to provide 
a homogeneous space charge region. However, for the 
realistic parameters = 10 15 cm~ 3 and a = 30 nm (so 
a <C l c ) the volume in which the potential drops to half 
its initial value contains approximately one doping atom. 
Discrete energy levels of such a doping atom cannot be 
resolved at room temperature. More importantly, the 
potential well due to an ionized single dopant will locally 
distort the barrier shape. This effect complicates the 
potential landscape, but it can only significantly increase 
the conduction of the diode, when the dopant resides 
close to the interface jlfj . 

In conclusion, we have shown by means of a simple 
electrostatic argument that the Schottky barrier thick- 
ness becomes a function of the diode size for small diodes 
(e.g. smaller than l c w 80 nm for TVa = 10 17 cm -3 ). Con- 
sequently, the contribution of tunneling to the total con- 
ductance is greatly enhanced in small diodes. This effect 
explains several experimental results (F^, 0], without the 
assumption of a reduced SBH. Moreover, small diodes 
show TV-curve shapes that qualitatively differ from those 
of conventional diodes. 
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